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We investigate the transport characteristics of Ge/Si-Core/Shell nanowire with Coulomb Blockade 
in presence of external magneto-electric fields from a theoretical basis. Using the effective Luttinger- 
Kohn Hamiltonian we calculate the valence band energy states of the nanowire and find out the state 
density at mili-Kelvin temperature. We explore the current-voltage and conductance-voltage char- 
acteristics. The observed transport characteristics are in close agreement with experimental data. 
We find these characteristics to be sensitive to the coupling of the quantum dot with the reservoirs. 
The quantum nature of the characteristics becomes less prominent with increasing temperature, as 
expected. 

PACS numbers: 73.23.Hk,73.63.Kv,73.63.-b 



I. INTRODUCTION 

Engineering spin splitting of charge carriers in 
semiconductor nanostructures such as quantum 
wires and dots has opened up realistic opportu- 
nities for the production, detection and manipu- 
lation of spin currents, allowing coherent trans- 
mission of inf orm ation beyond classical informa- 
tion processin^n^l jt is expected that spatial con- 
finement will have a significant effect on Zeeman 
splitting when bound-state quantization energies 
are comparable in nano-devices. Recent advances 
in fabrication technology have created opportuni- 
ties to investigate hole spin physics in semiconduc- 
tor nanowires for a range of different material^^. 
To achieve the proposal of Loss and DiVincenzcP, 
scientists thus have fabricated quantum dots in- 
tegrated with charge sensors namely as quantum 
point contaclP to initialize, manipulate and read- 
out electron spins in GaAs system J^l. However, the 
attempt suffers from difficulties in control and lack- 
ing of spin coherence due to the overlap of the elec- 
tronic wavefunctions in these systems with a large 
number of nuclear spins^. Though this intrinsic 
decoherence rate has been successfully reduced by 
spin-echo techniqu^, the complexity of the system 
remains as an issue under debate. 

As a result, the prospect of long coherence 
times in group IV materials due to the predom- 
inance of spin-zero nuclepsl has stimulated sev- 
eral proposalJii' and significant experimental ef- 
fort. Experimental progress includes realizations 
of DQDs in carbon nanotubed^ and SirF^^ as 
well as single dots in Si and Ge/Si nanowirei^l 
and Si/Ge 2DEG^. Moreover, the observation 
of long relaxation times due to negligible inter- 



action with nuclear spin bath makes the spin of 
a hole in a semiconductor quantum dot the best 
contender for solid-state qubit. Thus, Ge/Si- 
Core/Shell nanowire has become a central attrac- 
tion for both the experimental and theoretical in- 
vestigation since the valence band offset at core- 
shell interface provides a strong radial confinement 
of holes and due to the large sub-band spacing, this 
behaves as one-dimensional hole gas at ultra low 
temperatur^^. 

In this work, we explore the Coulomb Blockade 
effect in current-voltage and conductance-voltage 
characteristics of Ge/Si-Core/Shell nanowire 
Quantum Dot under perpendicular electric and 
magnetic fields. The Valence band states of 
the nanowire are theoretically analyzed us- 
ing Luttinger-Kohn HamiltoniarP^HH Then 
a self-consistent field approximation using 
non-equilibrium Greens Function formalism is 
performed to simulate density of states and from 
there we investigate the transport characteristics 
in mili-Kelvin temperatures. 



II. MODEL HAMILTONIAN AND 
SELF-CONSISTENT SOLUTION 

In our present analysis we use an effective 
ID Luttinger-Kohn Hamiltonian for a cylindrical 
nanowire with Ge core of radius, i? = 5 nm and Si 
shell of thickness t — 2 nm. This Hamiltonian is 
perturbative and corresponds to ground state and 
first excited state only. The electrochemical po- 
tentials of the materials are comparable to the en- 
ergy of these states and the applied potentials are 
in mili-volt order, hence the higher excited states 
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will not contribute significantly in current conduc- 
tion. For our calculation we assume the confine- 
ment to be an infinite cylindrical well, since the 
valence band offsets in Si/Ge Core/Shell nanowire 
is large. Now, the four eigenstates, \g±) and |e±), 
corresponding to ground and excited states for 
Fz — ±1/2 {Fz is quantum number such that, 
[Hlk + V, Fz] = and the operator Fz — Lz + Jz, 
Lz being orbital angular momentum along the wire 
axis) sX kz —0 will serve as the basis states for the 
Hamiltonian. Here, we have added the contribu- 
tion of electric field that causes Rashba Spin Or- 
bit InteractiorP^, the effect of external magnetic 
field which lifts the Kramer degeneracy and the 
spin splitting due to strain which is induced by 
the compression of Ge core by Si shell. A detailed 
derivation of the Hamiltonian for G e/Si core/shell 
nanowires is provided elsewher d^^ l ^^ ^ here we quote 
the matrix needed to calculate hole spectrum. 
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h^kf_ 0.73/12 
H22 = H44 = —^ + + <5(7) 

H12 = H34 = -I —kz 

moK 

Hi3 = ^lbBx{Xi + X2)] H24 = l^BBxiXi — X2) 
Hi4 = -~H23 = -i^isBxXs - ieEyUo 

with, ^ 10 T, Ey = 1 V/m. For Ge, 

Uo = 0.15R, X3 = 8.04i?. Xi and X2 are ma- 
terial dependent constants and 5{'j) rescales the 
energy splitting and depends on the relative shell 
thickness^^, 7 — t/R. The energy eigenvalues are: 
-4.675 meV, -1.45 meV, 5.075 meV, and 8.15 meV. 



III. DENSITY OF STATES 

We denote the creation and annihilation opera- 
tors of the ground state and the first excited state 
by c, c+ and d, respectively. The Schrodinger 
equation for these operators incorporating the 
Coulomb interaction between them becomes: 



in—c ■■ 
At 



Ud+dc 



(2) 



Here e is the energy of the energy level c, is the 
coupling between the level c of the channel and 



the r-th energy level of the reservoir, and is the 
energy of the r-th reservoir energy level. U is the 
coulomb potential due to interaction between the 
two energy levels of the channel. We write the 
uncoupled reservoir operators as Cr and define a 
source term S"^ := Cr in Fourier space for the 
channel operators. Then the Fourier space solution 
of ^ can be written as c = G{E)S'' where G{E) 
is the Green's functiorP^: 
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The density of states, denoted by D{E), is calcu- 
lated using self-consistent approximation using the 
Green's function ([s]). We find. 



D{E)^i{G-G* 



(5) 



The external magnetic field breaks the degeneracy 
of both the ground state and the first excited state 
and each of these states separates into two different 
states corresponding to spin up and down. 



IV. I-V CHARACTERISTICS 

We assume that the temperature is low enough 
for particle production due to thermal excitation 
to be ignored. Therefore amount of charge in the 
channel at any time is well defined. We also assume 
a constant capacitance as in the constant interac- 
tion modeP^. From the density of state ^ the 
current / is easily evaluated, 

I{Vg,Vsd) ^ -rf^, 

Kii +72) 

X / dED{E-U -VG){fs{E + VsD/1) 

J 00 

-fD{E-VsDm) (6) 

where fs and fu are Fermi functions of the source 
and drain. 

First we observe the variation of current with 
varying gate voltage (Vg) at constant source-drain 
voltage (Vsd)- In Fig- [l] (top) discrete cur- 
rent peak appears, meaning Coulomb blockade is 
lifted by changing the gate voltage. When the 
electrochemical potential of the dot in presence 
of one hole, is in the bias window, in our case 
(AiSiMu) — (0.2 meV, 0.3 meV), one hole tunnels 
into the dot from the source, so that the number of 
holes increases from to 1, causing the first peak at 
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-6.4 mV. This single hole tunneling continues and 
the electrochemical potential level corresponding 
to transport rises sequentially and the number of 
holes increases upto 4 for Vg — —3.25 mV, 3.2 mV 
and 6.5 mV. In the valleys between these peaks, 
the numbers of holes on the dot is fixed due to 
Coulomb Blockade. By tuning the gate voltage 
from one valley to the next one, the number of 
holes on the dot thus can be precisely controlled. 
The distance between the peaks corresponding to 
-E'add which in our case are found 3.15 meV for 
ground spin states and 3.3 meV for excited spin 
states respectively. 
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FIG. 1. (Top) Variation of current with respect to 
gate voltage. Here, Blue: Vsd ~ 0.5 mV, Red: Vsd ~ 
1.5 mV, Green: Vsd ~ 3 mV. (Bottom) Variation of 
current with respect to drain-to-source voltage. Here, 
dashed: Vg ~ mV, Green: Vg = —2.5 mV, Brown: 
Vg = -5 mV. 



cited spin states respectively. Moreover, increas- 
ing Vsd signifies transition involving excited state 
falls within the bias window and double hole tun- 
neling occurs which implies the number of holes 
alternates between 0, 1 and 2 leading to a double 
hole tunneling. According to Fig. [l] (top) , it can 
be concurred that the peak current increases from 
approximately 2.35nA to 3.18nA for the variation 
of source to drain voltage from 0.5 mV to 3 mV. 

In Fig. [T] (bottom) when the voltage is suffi- 
cient for a states of iV = 1 hole to be transferred 
to the dot to the reservoir, a sudden jump of the 
current will first be observed due to the junction 
with a high tunnel rate. Then the current will 
grow slowly, because of the low-rate junction com- 
ing into play, until states of iV = 2 hole appears 
on the dot. Thus, although the current through 
the system passes continuously, a particular volt- 
age dependent amount of holes will exist showing 
" Coulomb staircase" . 
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FIG. 2. Contour plot of I-Vg-Vsd- Here, transport 
is in z-a^xis, radius of Ge-core R — 5 nm, Rashba SOI 
Ey — 1 V/m, applied magnetic field, Bx ~ 10 T, tem- 
perature = 500 mK. 



We now look at the regime where Vsd is so 
high that multiple dot levels can participate in 
hole tunneling. A positive Vsd increases the dif- 
ference of electrochemical potentials of the source 
and the drain. Due to capacitive coupling the lev- 
els of the dot are also increased. In our case we 
assume, Cg = 3 aF, Cd = Cs = 10 aF. Thus, 
for Vsd = 1-5 mV, the hole states are separated 
further and Fig. [T] (top) shows that the number 
of holes increases from to I at -7.1 mV, 1 to 2 
at -3.85 mV, 2 to 3 at 2.5 mV and 3 to 4 at 5.9 
mV respectively. The distance between the peaks 
corresponding to i?add which are found to be 3.25 
meV for ground spin states and 3.4 meV for ex- 



In Fig. [2] we plot a contour of current in the 
Vq-Vsd plane as a representation of the I-V char- 
acteristics. Here, each time a hole is added to the 
quantum dot, a current peak appears as a diagonal 
line, a typical signature of a single quantum dotPsl. 



V. CONDUCTANCE CHARACTERISTICS 

Now from the expression of current (|6| we eval- 
uate the conductance, 

aiVa, Vsd) = J^—HVg, Vsd)- (7) 

OVsD 
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The conductance response of the device completely 
characterizes it and here we observe the contours 
of conductance in the Vg-Vsd plane. 

Conductance (nS) contour 
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FIG. 3. Contour plot of conductance- Vg-Vso. Here, 
transport is in z-axis, radius of Ge-core R = 5 nm, 
Rashba SOI Ey = 1 V/m, applied magnetic field, Bx — 
10 T, temperature = 500 mK. 

To understand Fig. [3j first the electrochemical 
potential ladder is defined by the gate voltage axis 
of the {VsD and Vq) plot. Now, by sweeping the 
gate voltage at low bias regime, hole tunneling oc- 
currs due to the transition \g^) o \g+). For all 
other gate voltages the dot is in Coulomb block- 
ade. The figure also depicts several sets of bright 
lines. These are the regions of maximum differen- 
tial conductance, i.e., the dot current changes in 
these sets of Vsd and Vq- The slopes of the V- 
shaped lines at 6 mV depend on the assumed ca- 
pacitance and thus defines the regions of Coulomb 
Blockade. The other solid lines indicate where the 
current changes due to the transitions involving 
excited states. From Fig. [3] it is evident that the 
separations of the ground and excited states for a 
single hole are 3 mV, 9.5 mV, and 12.75 mV. 

VI. EFFECTS OF TEMPERATURE 

In this section, we present a comparative study 
of the temperature dependence of hole transport. 
We simulate the self-consistent charge model in 
coulomb blockade regime in two different temper- 
atures, 500 mK and 5.0 K, and obtain the corre- 
sponding I-V curves and contours of conductance- 
voltage. With increasing temperature, it is evident 
that the transport should become more classical 
and hence quantum behavior of the system should 
be less prominent. 

In Fig. |4j one can see that the peaks of current 
with respect to gate voltage are smearing out with 
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FIG. 4. Variation of current with respect to gate volt- 
age for T = 500 mK (top) and 5.0 K (bottom). Here, 
Blue: Vsd = 0.5 mV, Red: Vsd = 1.5 mV, Green: 
Vsd = 3 mV. 
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FIG. 5. Variation of current with respect to drain- 
to-source voltage for T — 500 mK (top) and 5.0 K 
(bottom). Here, dashed: Vq = mV, Green: Vg = 
—2.5 mV, Brown: Va = —5 mV. 



increase in temperature. This is because energy 
levels become excited due to thermal fluctuations 
and thus discreteness of the curve tends to vanish. 
In Fig. [5j we have plotted current with respect 
to drain-to-source voltage. Due to thermal excita- 
tion, the staircase tends to be linear and becomes 
dispersed with increasing temperature. 

Finally, in Fig. [6) contours of current / and 
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FIG. 6. Contour plot of I-Vg-Vsd (top) and 
conductance- Vb-VsD (bottom) for T = 500 mK (left) 
and 5.0 K (right). Here, transport is in z-axis, radius of 
Ge-core R — 5 nm, Rashba SOI Ey = 1 V/m, applied 
magnetic field, Bx ~ 10 T. 



conductance a in the Vg-Vsd plane are plotted 
and they also lose their discrete nature with the 
rise in temperature from 500 mK to 5.0 K. All 
these results are physically consistent and can be 
explained by thermal excitation. 



VII. DISCUSSION 

In our analysis we have been able to re- 
produce most the prominent features of 
the Ge/Si -Core /S hell nanowire as found in 
experiment j^^ l ^^ l^ and the developed model is 
capable of describing most of the detailed features 
by taking into account the properties of the 
materials. The self-consistent field approximation 
using non-equilibrium Greens function formalism 
simulates realistic density of states and transport 
characteristics in mili-Kelvin temperatures. We 
observe the quantization of I-Vq and the staircase 
profile of I-VsD- Additionally, the contour of 
Conductance- Vg-^sd characterizes the device 
operating in single hole transition. We find these 
characteristics to be significantly sensitive to the 
coupling of the dot to the reservoir by means 
of capacitance effect. They are also sensitive 
to the operating temperature and start losing 
their quantum nature with rise in temperature. 
Hence, from a theoretical point of view, this work 
suggests that Ge/Si-Core/Shell nanowire is an 
attractive candidate for quantum information 
processing. However, several directions seem 
promising for further explorations. For this 
particular device structure, even cleaner wires 
would show evidence of asymmetric geometrj^SHI] 
and should be accounted for in theoretical model. 
Moreover, soft wall consideration along with wave 
function penetration can be added to achieve 
results with better accuracy. 
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